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Abstract 

While nanocatalysis is a very active field, there have been very few studies in the size/shape-dependent catalytic 
properties of transition metals from a thermodynamical approach. Transition metal nanoparticles are very attractive 
due their high surface to volume ratio and their high surface energy. In particular, in this paper we focus on the 
Pt-Pd catalyst which is an important system in catalysis. The melting temperature, melting enthalpy, and catalytic 
activation energy were found to decrease with size. The face centered cubic crystal structure of platinum and 
palladium has been considered in the model. The shape stability has been discussed. The phase diagram of 
different polyhedral shapes has been plotted and the surface segregation has been considered. The model predicts 
a nanoparticle core rich in Pt surrounded by a layer enriched in Pd. The Pd segregation at the surface strongly 
modifies the catalytic activation energy compared to the non-segregated nanoparticle. The predictions were 
compared with the available experimental data in the literature. 

PACS: 65.80-g; 82.60.Qr; 64.75.Jk 



Introduction 

Bimetallic nanoparticles exhibit unusual physicochemical 
properties different from those of the bulk material or 
their individual constituents [1,2]. They are very used in 
catalysis, fuel cells, and hydrogen storage. These unusual 
properties are determined by their size, shape, and com- 
position. When considering metallic catalysts, platinum 
is a standard material but this material is most expen- 
sive than gold [3]. Therefore, to reduce the amount of 
platinum and then the cost of the application, one possi- 
ble way is to use an alloy of platinum with another 
metal. In the present study, the chosen alloy is the bin- 
ary Pt-Pd system [4] that we propose to theoretically 
study from a thermodynamic approach [5,6], as well as 
its pure components. It has been shown previously [5,6] 
that thermodynamics may provide useful insights in 
nanotechnology where the size of the considered nano- 
particles is higher than approximately 4 nm. Within this 
approach, the size and shape effects on the melting tem- 
perature, melting enthalpy, phase diagram, and catalytic 
activation energy of this system are investigated. 
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As face-centered cubic (fee) metals, Pt and Pd can 
exhibit a variety of geometrical shapes. Therefore, to 
address the shape effect on the materials properties of 
these metals at the nanoscale [7,8], the following shapes 
have been considered: sphere, tetrahedron, cube, octahe- 
dron, decahedron, dodecahedron, truncated octahedron, 
cuboctahedron, and icosahedron. 

Size-dependent melting properties of Pt and Pd 

At the nanoscale, the melting temperature T m and melt- 
ing enthalpy AH m , for free-standing nanostructures can 
be expressed as function of their bulk corresponding 
property, the size of the structure and one shape para- 
meter [9]. 



^m/^ra,oo - 1 ^shape/^A 
AH m / AH moQ = T m J Tm,ooi 



(i) 



(2) 



where the shape parameter, a shape , is defined as a shape 
= AD(y S 'Yi)/(VAH m} J; D being the size of the structure 
(i.e. for a sphere, D is the diameter), A (meter squared) 
and V (cubic meter) are the surface area and volume of 
the nanostructure, respectively. AH m>00 is the bulk melt- 
ing enthalpy (Joule per cubic meter), whereas // and y s 
are the surface energy in the liquid and solid phases 
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(Joule per square meter), respectively. // and y s are con- 
sidered size independent. This is justified by the fact 
that the size effect on the surface energies is less than 
4% for sizes higher than 4 nm [10,11]. Indeed, below 
this size, edges, and corners of the structures begin to 
play a significant role in the surface energy [12]. 

The size-dependent melting temperatures of platinum 
and palladium are plotted in Figures 1 and 2 respec- 
tively. The materials properties of the considered mate- 
rials are indicated in Table 1. The melting properties for 
the sphere have been calculated using for the solid sur- 
face energy the mean value of experimental data [13]. 
For the other polyhedra shapes, we have considered the 
fee crystal structure of the metals and the respective 
solid surface energy for each face [14]. Tables 2 and 3 
indicate the parameters used for the calculation of the 
melting properties. Experimentally, the melting of 
agglomerated Pt nanocrystals (tetrahedrons and cubes) 
with an average size around approximately 8 nm starts 
at approximately 900 K [15] in relative good agreement 
with our theoretical predictions. Molecular dynamics 
simulations [16] have calculated the size effect on the 
melting temperature of Pd and found a spher e = 0.95 nm 
while our theory predicts 1.68 nm. 

Discussion 

At the nanoscale, the shape which exhibits the highest 
melting temperature is the one which minimizes the 
most the Gibbs' free energy (G = H - TS); and is then 
the favored one. From Figures 1 and 2, the four most- 
stable shapes among the ones considered are the dode- 
cahedron, truncated octahedron, icosahedron, and the 
cuboctahedron. Experimentally, truncated octahedron 
and cuboctahedron are observed for platinum nanopar- 
ticles [8] whereas icosahedron, decahedron, truncated 
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Figure 1 Size-dependent melting temperature of platinum 
versus the size for different shapes. 
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Figure 2 Size-dependent melting temperature of palladium 
versus the size for different shapes. 



octahedron and cuboctahedron are observed for palla- 
dium nanoparticles [8]. Therefore, our predictions are in 
relative good agreement with the observations for palla- 
dium and platinum except that dodecahedron and icosa- 
hedron are not observed for platinum. Other theoretical 
calculations confirmed that the dodecahedron is a stable 
shape for palladium [17]. More generally, according to 
Yacaman et al [8], the most often observed shapes at 
the nanoscale are the cuboctahedron, icosahedron, and 
the decahedron. 

Furthermore, care has to be taken when we compare 
theoretical results with experimental ones due those 
materials properties depend on the synthesis process 
[18,19]. And then predicted properties from thermody- 
namics may differ from the experimentally observed if 
the synthesis process is not running under thermodyna- 
mical equilibrium. Moreover, thermal fluctuations are 
often observed in nanoparticles [20] meaning that the 
shape stability is much more complicated than just a 
minimisation of the A/V ratio with faces exhibiting the 
lowest surface energy. 

Nano-phase diagram of Pt-Pd 

According to the Hume-Rothery's rules, platinum and 
palladium forms an ideal solution [21]. In this case, con- 
sidering no surface segregation, the liquidus and solidus 

Table 1 Materials properties of platinum and palladium 



Materials properties 


Platinum 


Palladium 


W (K) [40] 


2,041.5 


1,828 


AH m/OQ (kJ/mol) [40] 


22 


17 


AH sub , M (kJ/mol) [41] 


565 


377 


7/ (J/m 2 )[40] 


1.866 


1.470 


7s (Vm 2 ) [13] 


2.482 


2.027 
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Table 2 Solid surface energies for platinum and 
palladium materials [13]. 



Faces 


Platinum 


Palladium 


Y (111) (Vm 2 ) 


2.299 


1.920 


y 5 (100) (J/m 2 ) 


2.734 


2.326 


^ (110) (J/m 2 ) 


2.819 


2.225 



curves of bulk and nanostructures are calculated from 
the following equations [22-24]: 



(■^solidus 
■^liquidus 
1 ^solidus 



AHi 



\ m 



V 1 ^liquidus / \ T m ) 



(3) 



where # so iidus (#Hquidus) is the composition in the solid 
(liquid) phase at a given I 7 , respectively. T l m is the size- 
dependent melting temperature of the element i. AH l m is 
the size-dependent melting enthalpy of the element i. 

The phase diagram of the Pt-Pd alloy is plotted in Figure 
3. We note that the lens shape of the phase diagram is 
conserved at the nanoscale; however, the lens width 
increases for the shapes characterized by a small melting 
enthalpy and melting temperature, i.e., exhibiting a strong 
shape effect. Moreover, the melting temperature increases 
with the concentration of Pt in agreement with Ref. [25]. 

In order to predict nanomaterials properties more 
accurately, we are considering a possible surface segre- 
gation which is known as the surface enrichment of one 
component of a binary alloy. At the nanoscale, surface 
segregation leads to a new atomic species repartition 
between the core and the surface. According to Wil- 
liams and Nason [26], the surface composition of the 
liquid and solid phase are given by: 



^surface 
solidus 



surface 
liquidus 



(^solidus / 1 ^solidus) # 



-(AH suh z lv )/( Zl kT) 



1 + (Xsolidus/1 " Solidus) erW^ltoKT) 
(Xliquidus/l ~ ^liquidus) e~( AH ^/^ 
' 1 + (Xliquidus/l " liquidus) e-( AH ^/^ 



,(4) 



Table 3 Number of (hkl) faces for each shape. 



Shape 


Number of 
(111) faces 


Number of 
(100) faces 


Number of 
(110) faces 


Tetrahedron 


4 


0 


0 


Cube 


0 


6 


0 


Octahedron 


8 


0 


0 


Decahedron 


10 


0 


0 


Dodecahedron 


12 


0 


0 


Truncated 
octahedron 


8 


6 


0 


Cuboctahedron 


8 


6 


0 


Icosahedron 


20 


0 


0 




Figure 3 Phase diagram of the Pt-Pd system for different 
shapes. Different shapes at a size equal to 4 nm and at the bulk 
scale. The solid lines indicate the liquidus curves whereas the 
dashed lines indicate the solidus ones. 



where z x is the first nearest neighbor atoms; z lv is the 
number of first nearest atoms above the same plane 
(vertical direction). In the case of face-centered cubic 
(fee) crystal structure of Pt and Pd materials, we have Z\ 
- 12, Zi v - 4 for (100) faces and three for (111) faces. 
A// vap is the difference between the bulk vaporization 



enthalpies of the 



AH V: 



AH" 



pure 
is the 



elements, 
difference 



two 

A Z_7 

I vap - AA1I y,oo AA1I y,oo' ^ ''sub 

between the bulk sublimation enthalpies of the two pure 
elements, AH su b = AH^ — AHf 0G . Element A is chosen 
to be the one with the highest sublimation and vaporiza- 
tion enthalpies. If the two components are identical, 
A// su b = 0 and A// vap = 0, there is no segregation and 
we retrieve Equation 3. # so iidus an d #ii qu idus are obtained 
from solving Equation 3. Assuming an ideal solution, 
only the first surface layer will be different from the 
core composition. 

Considering the surface segregation in the Pt-Pd sys- 
tem, we can see in Figure 4 that the lens shape of the 
surface liquidus/solidus curves is deformed compared to 
the core. At a given temperature, the liquidus and soli- 
dus curves of the surface are enriched in Pd compared 
to the core; meaning that the surface is depleted of Pt 
(the higher bond energy element) which is in agreement 
with experimental observations [27-29] and other theore- 
tical calculations [29-31]. This is due to the fact that Pd 
has a lower solid surface energy, a lower cohesive energy 
compared to Pt and also because diffusion is enhanced 
at the nanoscale [32]. 

Size-dependent catalytic activation energy of Pt-Pd 

The catalytic activation energy is the energy quantity 
that must be overcome in order for a chemical reaction 
to occur in presence of a catalyst. The low the catalytic 
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Figure 4 Phase diagram of the Pt-Pd system considering the 
surface segregation effect. Surface segregation effect at a size 
equal to 4 nm for a spherical nanoparticle. 
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Figure 5 Composition dependency of the catalytic activation 
energy for a spherical nanoparticle of Pt-Pd. Nanoparticle of Pt- 
Pd with a size equal to 4 nm. 



activation energy is, the most active the catalyst is. It is 
thus an important kinetic parameter linked to the che- 
mical activity. Indeed, the catalytic activation energy is a 
linear function of the work function [33-35]. For pure 
materials, the catalytic activity depends on the fraction 
of surface atoms on corners and edges while for binary 
compounds it depends also on the surface segregation. 
Recently, it has been showed by Lu and Meng in Ref. 
[36] that the size-dependent catalytic activation energy, 
£ ca could be obtained from the following relation: 

^ca / ^ca,oo = T m j T moo (5) 

Therefore, it means that the size-dependent catalytic 
activation energy decreases with size. 

To compare with experimental results, the ratio of the 
catalytic activation energies between tetrahedral (D = 
4.8 nm) and spherical (D = 4.9 nm) pure platinum 
nanoparticles has been determined around 0.66 in excel- 
lent agreement with the experimental value of 0.62 ± 
0.06 announced by Narayanan and El-Sayed [37-39]. 
Moreover, the ratio of the catalytic activation energies 
between cubic (D = 7.1 nm) and spherical (D = 4.9 nm) 
pure platinum nanoparticles is around 1.01 in relative 
good agreement with the experimental value of 1.17 ± 
0.12 [37-39]. 

From the size-dependent Pt-Pd phase diagram, the 
melting temperature of the alloy can be deduced. Equa- 
tion 6 describes the melting temperature of the bulk Pt- 
Pd while Equations 7 and 8 describe the nanoscaled 
melting temperature of a non-segregated and segregated 
spherical nanoparticle (with a diameter equal to 4 nm), 
respectively. 



Tsolidus (Bulk) = 1828 + 23 6x - 22x 2 , (6) 
T so iidus_core (D = 4 nm) = 1019 + 258x - llx 1 , (7) 

T so iidus_surface (D = 4 nm) = 1264- 111 exp (-x/0.01 6) -58 exp (-x/0.0020) -73 exp (-x/0. 1043) , (8) 

where x represents the alloy composition. For a sphe- 
rical Pt-Pd nanoparticle with a diameter equal to 4 nm, 
by combining Equations 5-8, £ ca seems to evolve quad- 
ratically with the composition when the segregation is 
not considered; which is not the case when the segrega- 
tion is considered (Figure 5). For the segregated Pt-Pd 
nanoparticle, a maximum in the catalytic activation 
energy is reached around 16% of Pt composition. 

Conclusions 

In conclusion, it has been shown that thermodynamics 
can still provide useful insights in nanoscience and 
more specifically in catalysis. The future development of 
catalysts and fuel cells is dependent upon our ability to 
control the size, shape, and surface chemistry of indivi- 
dual nanoparticles. Future theoretical work will have to 
consider the environment in which the particles are 
synthesized as well as the preparation method because 
these parameters can have a great influence on the 
shape stability and on the catalytic properties. 
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